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ABSTRACT 

We use Monte Carlo simulations to estimate the number of extrasolar planets that are di- 
rectly detectable in the solar-neighborhood using current and forthcoming high-contrast imaging 
instruments. Our calculations take into consideration the important factors that govern the 
likelihood for imaging a planet, including the statistical properties of stars in the solar neigh- 
borhood, correlations between star and planet properties, observational effects, and selection 
criterion. We consider several different ground-based surveys, both biased and unbiased, and 
express the resulting planet yields as a function of stellar mass. Selecting targets based on their 
youth and visual brightness, we find that strong correlations between star mass and planet prop- 
erties are required to reproduce high-contrast imaging results to date. Using the most recent 
empirical findings for the occurrence rate of gas-giant planets from RV surveys, our simulations 
indicate that naive extrapolation of the Doppler planet population to separations accessible to 
high-contrast instruments provides an excellent agreement between simulations and observations 
using present-day contrast levels. In addition to being intrinsically young and sufficiently bright 
to serve as their own beacon for adaptive optics correction, A-stars have a high planet occurrence 
rate and propensity to form massive planets in wide orbits, making them ideal targets. The same 
effects responsible for creating a multitude of detectable planets around massive stars conspire to 
reduce the number orbiting low-mass stars. However, in the case of a young stellar cluster, where 
targets are approximately the same age and situated at roughly the same distance, MK-stars 
can easily dominate the number of detections because of an observational bias related to small 
number statistics. The degree to which low-mass stars produce the most planet detections in 
this special case depends upon whether multiple formation mechanisms are at work. Relaxing 
our assumption that planets with large semimajor axes resemble the RV sample, our calculations 
suggest that the companions found in ultra-wide orbits around late-type stars in clusters are 
consistent with a formation channel distinct from that of RV planets. 

Subject headings: Stars: planetary systems - Methods: statistical, numerical 



1. Introduction 

Precision RV measurements have led to a statis- 
tically meaningful sample of extrasolar planet dis- 
coveries (http://exopl anets.org| that have shaped 
much of our understanding of planet formation 
and evolution. Presently, the RV-detected plan- 
ets found to orbit most distantly from their host 
star, 55 Cnc d a nd HP 190984 b, have semimajor 
axes of w 6 AU ( Fischer et al.ll2008t ISantos et al 



20101) . As the time baselines of RV surveys 
grow and more planets are found at incremen- 
tally longer periods, a complementary technique, 
high-contra st imaging, shows promise to suppl 
detections dMarois et all 120081: | K alas et ail 1200 



Lagrange et al. 2010: Marois et al. 20101) and 



et al.l 
20101) 



place tight constrain ts on the planet popul ation 
at large separations (jNielsen and Close 201dl ) by 
working from the opposing direction - outside-in. 
Of the planets detected directly, the two that orbit 
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closest to their host star, j3 Pictoris b and HR 8799 
e, have semi ma.jor axes of ~ 12 AU and ~ 15 AU, 
respe ctively (|Lagrange et al. 20101 : Marois et al. 
20101) . Given th e current precision s achieved by 



RV instruments (lMarcv et al.l l2008t iMavor et al 



2009) and rapid development of coronagraphy 
and adaptive optics technology ( Absil and Mawet 
20101 : lOppenheimer and Hinklevll2009l ). it is clear 
that this gap in parameter space will continue to 
narrow and gradually fill. 

While the prospects for continued discoveries 
are encouraging, it is interesting to consider the 
origin of planets detected thus far. Both of the 
aforementioned groups - planets orbiting close 
to their host star and those orbiting far away 
- have seem ingly formed by distinct mechanisms 
( Bolev 20091 ). One currently favored explanation 
is that core-accretion and other processes related 
to gravitational disk instabilities a re responsible 
for building each group r espectively dPollack et al 
1996; llda and Lirj|2004 lBosslll997l: iDurisen et al 



20071 ). This assertion is nominally based on semi- 
analytic calculations that indicate conditions are 
unfavorable for planets to form by core-accretion 
in situ at separations exceeding s=s 35 AU, because 
of a prohibitively slow aggregation of planetesi- 
mals (jDodson- Robinson et al.ll2009t iKratter et al. 
20101 ). Taken at face- value, this result implies that 



either: some of the HR 8799 planets have formed 
through a different channel than their brethren; 
the regions over which core-accretion and grav- 
itational instability operate can overlap and the 
two may possibly conspire; or that the resonant 
membe rs of a multi-plane t system can migrate in 



unison 



Crida et al.1 12009). The debate regarding 



how and where Jovian plane t s typically form is 
on-going dMarois et al. 2010 ; Currie et al. 2011 ; 
Ireland et all l201ll; iBossI l201ll) . These scenarios 



can be tested by studying the planet population 
at intermediate and wide separations. 

In addition to probing an unexplored parameter 
space, high-contrast imaging also provides access 
to photons arriving directly from the planet it- 
self and thus enables characterization studies at 
a le vel of detail comparable to transiting plan- 
ets ( Seager and Deming 2010h . Efforts to under- 
stand planetary atmospheres are currently un- 
derway, but theoretical models are currently in 
a primitive state, having only a handful of ob- 
jects with which to study and suffering from 



model i nput paramete r degen e racies dBowler et al 
2010bl; IJanson et all l201ll : ICurrie et al.l l201ll 
Madhusudhan et al.l 201 il k To further inform our 
theoretical conceptions of planet formation, evo- 
lution, and atmospheric physics, it is necessary to 
obtain a larger sample of companions that reside 
in orbits beyond re 10 AU. 

Knowing where to search can help as it may 
be possible to bias the target selection strate- 
gies of planned imaging surveys towards "planet- 
enriched" stars. For example, several RV pro- 
grams have conducted biased surveys by specif- 
ically targeting metal-rich stars, leading to a 
rapid increase i n the number of close-in planets 



rapic 

(lLaughlinll2000l;lFischer et al.ll2005t Ida Silva et al 
2006). If the number density of RV planets is 



growing with log orbital period, then a modest 
extrapolation to separations accessible to high- 
contrast imaging instruments may provide a rea- 
sonable starting point for maximizing the yield of 
future surveys. 

Recent studies of the Doppler RV population 
show that the planet occurrence rate scales not 
only with metallicity, but also with stellar mass. 
Johnson et al.l (|2010al ) have shown that the frac- 
tion of stars with planets can be parameterized 
as a function of stellar mass, M*, and metallicity, 
[Fe/H], according to: 

/(M*, [Fe/H]) = fc(M,/M riO^ Fe / H ], (1) 

where k = 0.07 ± 0.01, a = 1.0 ± 0.3, and 
[3 = 1.2 ± 0.2. This result is the first to prop- 
erly isolate the effects of stellar mass and metal- 
licity on the known planet population and could 
have important implications for current and fu- 
ture direc t-imaging program s, such as the VLT 
SPHERE dPohlen et al.ll2006|) Gem ini Planet Im- 
ager (GPI) (iMacintosh et al ll2006h . Project 1640 
at Palo mar ( Hinklcv et al. I~ l2011bh . and Subaru 
SEEDS (|Tamural l2009). These instruments cover 
a relatively large range of declinations, and com- 
bined, they will observe more than 500 stars in the 
solar neighborhood over the next several years. It 
would be useful to identify which kinds of stars are 
most promising to target. 

Simulations can be used to estimate the de- 
gree to which certain spectral-types are likely to 
host imageable planets. In addition to planet oc- 
currence rates, effects that must be taken into 
consideration include: the stellar mass function, 
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distribution of local stellar ages and metallicities, 
planet orbital properties, correlations between 
stellar properties and planet properties (such as 
mass and semimajor axis), near-infrared contrast 
ratio between the star and planet, and, in practice, 
wavelength range used by the science instrument 
and AO system wave-front sensor. With exception 
of the planet semimajor axis distribution, many 
of these factors are now either known empirically 
or based on well-constrained observational or the- 
oretical grounds: those related to the number 
statistics of stars and RV exoplanets in the solar 
neighborhood are backed by a wealth of high res- 
olution spectroscopic measurements; observation- 
related trade-offs are appreciated from previous 
experience with coronagraphic and AO hardware; 
and, although exoplanet atmospheric models are 
uncertain at young ages, their systematic over- 
estimate or under-estimate of the intrinsic bright- 
ness of planets in a given band may still be used 
to calculate direct detection numbers and rates in 
a relative sense. 

In th is paper, we use the results from ljohnson et al 
(|2010al ) as input to Monte Carlo simulations to 
estimate the number of planets that are directly 
detectable from the ground at near-infrared wave- 
lengths. We consider high-contrast imaging obser- 
vations that survey target stars using several dif- 
ferent representative samples. The results, which 
are based on the current understanding of planet 
formation as informed by observations, are tabu- 
lated as a function of stellar mass. They may be 
used to justify and guide target selection for di- 
rect imaging programs, in order to maximize their 
return, and to help understand the theoretical rel- 
evance of planet discoveries made in the future at 
intermediate and large orbital separations. 

2. Description of Monte Carlo Simulations 

To calculate the number of planets detectable 
by future surveys, we begin by synthesizing a 
volume- limited sample of stars within 50 pc of the 
Sun f £|2.ip . This distance is large enough to jus- 
tify the use of statistically significant, empirical 
results for the properties of stars within the solar 
neighborhood, yet close enough to yield a reason- 
able chance for detecting a planet directly with a 
large-aperture telescope in the near-infrared. It 
also corresponds to a space-volume wherein all 



A, F, G stars are bright enough at visible wave- 
lengths to serve as their own natural guide-star. 
We then populate the stars with planets, accord- 
ing to Eqn. [TJ and other parameters f H2.2[) . and 
simulate observations with a coronagraph and an 
"extreme" AO system fffO]). 

Several different large-scale surveys are consid- 
ered: a full volume-limited survey where each star 
within 50 pc is observed f H3.1j) : a survey in which 
stars within this volume are selected by age, to 
have relatively bright planets, and stellar visual 
magnitude, for effective AO correction ( §3.2j) ; a 
survey where the selection criterion is modified 
to incrementally include more distant stars (as in 
the case of a magnitude-limited sample); and ob- 
servations of a nearby association that resembles 
Tucana-Horologium ( §3.5|) . We also perform cal- 
culations using contrast levels comparable to those 
achieved by current instruments. 

2.1. Simulated Stars 

We consider pre-main-sequence, main-sequence, 
and subgiant stars with mass 0.4 < M*/Mq < 2.6 
as potential targets. Using a numb e r density 
of 0. 11 stars per pc 3 ( Reid et al. 20021 : lLatvshev 
19781 ) , we simula t e a to tal of 10,081 potential tar- 
gets (|Grav et all (|2003l ): NStED). This figure as- 
sumes a multiplicity factor of 0.33 for the fraction 
of systems that are binary (Ra ghavan et al.ll2010f ) 
and that ^10% of those systems have physical 
separations amenable to high-contrast imaging, 
includi ng the individual componen ts of wide bi- 
naries ( Holman and Wiegertl 1 19991 ). Each target 
is assigned a mass, age, and metallicity that con- 
forms to the known bulk statistical properties of 
stars in the solar neighborhood^ 

1. Masses in the range 0.4 < M*/M < 
2.6 (M2-A1 V-IV) are drawn from either 
a present-d ay-mass-function (PDMF) for 
field stars ( Reid et all |2002| ), or an ini- 
tial mass function (IMF) for stellar clus- 



Ideally, one would acquire stellar properties for actual stars 
in the solar neighborhood and base a study on their indi- 
vidual properties. However, any attempt to do so would 
draw upon a heterogeneous collection of measurements that 
would be incomplete for the volume we wish to consider. 
To ensure that our estimates of relative detectability as a 
function of stellar mass are not dominated by systcmatics 
or predisposition to certain stellar-types, our study draws 
upon a synthetic sample. 
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ters ( Miller and Scald Il979h . depending on 
the simulation. The probability density 
for masses is given by a broken power-law , 
dN*/dM* oc Ml. The IReid et all (l2002h 
PDMF has indices 7 = -1.35 for M* < 
l-O Mw and t = -5 2 for M. > 1.0M©. 
The iMiller and Scalol l| 19791 ) IMF has in- 
dices 7 = -1.25 for M* < 1.0M©, 7 = -2.0 
for 1.0 < M*/M < 2.0, and 7 = -2.3 for 
M* > 2.0Mq. These numbers agree well 
with NStED database queries. 

2. Ages in the range, 4 Myr < t ago < 13.7 Gyr, 
are drawn from a histogram probability dis- 
tribution that follows a volume-limited sam- 
ple of well-chara cterized FG stars withi n 
40 pc based on INordstrom et al. (2004). 
Fig. [T] shows the input used for our code. 
We assume that the age distribution of all 
spectral-types follows the same probability 
distribution. Maximum age constraints are 
applied to ensure consistency with main- 
sequence or subgiant evolution for a given 
star mass. 

3. Metallicities in the range, —1.0 < [Fe/H] < 
0.55, are drawn from a sliding-Gaussian 
probability distribution with dispersion, 
0[Fe/H] = 0.18 dex, and center, [Fe/H] , 
that changes as a function of age, i ago , ac- 
cording to [Fc/H] = 0.177 - 0.040 i ago , in 
order to si mulate the a g e-met allicity relation 
derived by IReid et al.1 (|2007h . owing to the 
fact that young stars tend to be metal-rich. 

4. Distances from the Sun, s, are drawn from a 
dN* = Atts 2 ds relation. 

We construct «10 3 ensembles (10,081 stars per 
ensemble) and average the number of planets de- 
tected in 0.1 Mq wide stellar mass bins to converge 
to results accurate to 0.1 planets per bin. 

2.2. Simulated Planets 

Gas-giant planets are assigned to targets based 
on the star's mass and metallicity according to 
Equ. [U We assume that 28% of host stars have 
multiple planets, with eithe r 2,3,4 companion s 
each based on random chance ( Wright et al.l 2009). 
For each planet, we generate orbits using Monte 
Carlo techniques that take into account projection 
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Fig. 1. — Stellar age distribution used for Monte 
Carlo simulations. Values are drawn for 0.4 < 
Mq < 2.6M & stars using constraints to ensure 
consistency with main-sequence and subgiant evo- 
lution. Figure adapted from the results of a kinc- 
matically un biased 40 pc volume-lim ited survey of 
FG stars by INordstrom et al.l (|2004l ). 



effects (ICarson et al l J2006t iNielsen et al.l l2008t 
iThalmann et al.ll2009t iBeichman et alj|2010h : 



Orbits are first constructed in a two di- 
mensional plane from a given semimajor 
axis and eccentricity. The true anomaly is 
drawn from a probability distribution that 
is weighted by the amount of time spent 
by the planet at a given location in the el- 
lipse. Then, the argument of pericenter and 
ascending node are chosen randomly and 
an inclination, i, is drawn from a proba- 
bility distribution, dn p /di oc sin i, where 
i = represents face-on orbits. Finally, 
the orbit is oriented in three dimensions rel- 
ative to a distant observer using the ap- 
propriate rotational transfor mation matrices 
(jMurrav and Dermott]l2000l ). The angle be- 
tween the star and planet on the sky is cal- 
culated from the projected separation and 
stellar distance. 



Semimajor axes in the range a < a max are 
selected considering two different distribu- 
tions: dn p /da cx a~ 61 and dn p /da oc a . 
The first case is representative of a (con- 
ventional) planet population that has ex- 
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perienced negligible outward migration or 
planet-planet scattering, wit h an index given 



by th e relationship from iCumming et al 



(|2008l ). while the second is more applicable 
to planets that h ave either formed in-situ at 



large separation (|Bolevll2009t iKratter et al 



20inh or are dynamicall y inter-active with 
their disk or one another ([Scharf and Menou 
20091 : IVeras et all l2009h . We note that 
the exact form of the semimajor axis dis- 
tribution affects the absolute number of 
planets detected and less-so the relative 
number. Further, previous authors have 
found only a weak dependence of their 
Monte Carlo simulation results (e.g., for 
the occurrence rate of brown-dwarfs from 
observations) on t he assumed sem i major 
axis distribution (iLafreniere et al.l l2007t 



Metchev and Hillenbran d 2009; Beichman et al 



20101) 



3. The maximal extent of the semimajor axis 
distribution is scaled as a function of stel- 
lar mass, a max = (M,/2.5M Q ) 4 / 9 a, to ac- 
count for the theoretically motivated con- 
cept that the width of the zone where plan- 
ets may form increases wi th stellar mass 
( Kennedy and Kenvonll2008 ) . We later relax 
this assumption in i )3.5l to isolate its effect on 
the number of planet detections. Two rep- 
resentative cases are considered: a = 35, 120 
AU. The first corresponds to a distance be- 
yond which the core-accretion process is ex- 
pected to b e come inefficient (see however 
ICurrie et al. I (120111) 1 and the second corre- 
sponds app roximately to the location of Fo- 
malhaut b (jKalas et alJl2008h . 



4. The Johnson et al. ( 2010al ) study analyzed 
planets with semimajor axes a < 2.5 AU. 
The constant, k = 0.07, from Equ. [T] there- 
fore represents a lower-limit to the planet 
occurrence rate. To calculate the over- 
all occurrence rate, k e , we extrapolate the 
drip/da oc a~ 61 semimajor axis distribution 
to a = 35, 120 AU (corresponding to 23.3, 
79.9 AU for M* = 1.0M© stars), finding 
k e = 0.17,0.27 for each case respectively. 

5. We account for the observed paucity of plan- 
ets within 0.5 AU around A-stars by for- 
bidding semimajor axes in this region for 



stars with mass M* > 1.5M ^ (jBowler et al 
2010aH Johnson et alJTioiObh . 



6. Eccentricities are drawn from a probability 
density distribution, dn p /de = 1 — e. A lin- 
ear function that falls to zero when e = 1 
well-replicates the number statistics seen in 
data for exoplan ets with orbita l periods ex- 
ceeding 10 days (|johnsonll2009h . 



7. 



Planet masses in the range, 0.5 < m p /Mj < 
15, are drawn from a power-law distribu- 
tion, dn^/dm p oc mj, where 7 = —1.4 
(Johnson 2009). We assume that the in- 
dex of the true planet mass distribution 
equals that of the observed m p sin i dis- 
tribution. Statistical analyses suggest that 
this assumption is reasonable for isotropic 
orbital inclinat ions in the frame of the plan- 
etary system (jJorissen et al.l 120011 ). and in- 
direct constraints on the sin i distribution, 
via combined measurements of stellar rota- 
tional velocities, periods, and radii, confirm 
this result for well-aligned planetary sys- 
tems, i.e., those with small spin-orbit angles 
(jWatson et alj|201dh . 



Finally, we consider cases where A-stars 
have a propensity to form more massive 
planets than solar-type stars, by setting 
74 = —1.0, where is the power-law 
index as applied to planet masses around 
M* > 1.5M(D stars (jLovis and Mavorl [2007 : 



Johnson et al.l 2008). 



Table [T] summarizes the number-statistics used for 
both stellar and planetary properties. 

Planet masses are converted to absolute mag- 
nitudes in_e^£li_ofthe_^_= J, H, K bands us- 
ing th^jBarafe^etinjI^P^ "hot-start" (cond03) 
and lFortnev et al.l fljjjp ^ore-accretion" (fort08) 
evolutionary models. These models make dif- 
ferent predictions for the int rinsic bright n ess o f 
planets at young ages (see Marlev et al. ( 20071 ) 
for a discussion). We consider ages from 4 Myr 
to 13.7 Gyrs by interpolating between the values 
shown in the tables from the above references. To 
cover this entire range, we also extrapolate each 
model in time u sing the slopes for co oling curves 
from Fig. 1 of iFortnev et al which are 

roughly linear in magnitude vs. log(£ agc ) space, 
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Parameter 


Description 


Probability Distribution 


Reference 


Stars 


mass 


broken power-law 


7pdmf — 1.35, 5.2 


Reid et al. (2002) 






7lMF = -1.25,-2.0,-2.3 


Miller and Scab (1979) 


age 


histogram 


piece- wise, empirical 


Nordstrom et al. (2004) 


metallicity 


sliding Gaussian 


c [F c/H] = 0.18 dex 


Reid et al. (2007 s ) 


Planets 


mass 


power-law 


7mkgf = -1.4, 7a = -1-4, -1.0 


Johnson (2009) 


semimajor axis 


maximal extent 


a = 35, 120 AU 


Dodson-Robinson et al. (2009) 




power-law 


a max = (M»/2.5M©) 4 / 9 6 


Kennedy and Kenvon (2008) 




power-law 


7 = 0,-0.61 


Cummine et al. (2008) 


eccentricity 


linear 


dnp/de — 1 — e 


Johnson (2009) 



Table 1: Input parameters to Monte Carlo simulations. Power-law distribution indices are given by 7 for 
the equation dN/dx oc x 1 , where dN/dx is the number density as a function of the variable x. Stellar 
PDMF indices correspond to the mass range: M*/M < 1.0 and M„/M & > 1.0 respectively. Stellar IMF 
indices correspond to the mass range: M*/M© < 1.0, 1.0 < M*/M© < 2.0, and M* > 2.0M© respectively. 
The stellar age distribution is a histogram based on a volume-limited sample of well-characterized stars 
within 40 pes. The metallicity distribution is a function of age and has dispersion <7 [Fo / H ] = 0.18 dex. The 
maximal extent of planet semimajor axes, a = a, is scaled according to stellar mass. The planet eccentricity 
distribution is applicable to long-period planets discovered by the RV method. 



and s lopes for cooling curves from iBaraffe et al 
(|2003l ) , which are more appropriate for older plan- 
ets with low effective temperature s. We also ex- 
trapolate the iFortnev et al. ( 2008 ) core- accretion 
model to 0.5Mj and 15Mj based on the values 
provided in the range 1 < m p /Mj < 10. We 
do not consider the contribution to planet bright- 
ness from reflected starlight since thermal emission 
dominates at near-infrared wavelengths for young 
bodies with separations exceeding several AU, par- 
ticularly those that orbit low-mass stars. Stellar 
masses are conve rted into abso l ute m agnitudes in 
each band using Girardi et al.l (|2002[ ) isochrones. 
The magnitudes of the star and planet are then 
differenced to calculate contrast ratios. 

While the accuracy of planet evolutionary mod- 
els is a large uncertainty in the calculations pre- 
sented, we note that the choice of model is of sec- 
ondary importance because we are concerned pri- 
marily with the relative behavior in the number 
of planets detected as a function of stellar mass. 
In fj3l we show that the results from each model 
are qualitatively similar. We also note that the 
cond03 and fort08 models are representative of 
the extrema in brightness predictions at young 
ages and that the thermal evolution of a plane- 
tary atmosphere is likely to fall in between the 



cooling curves from "hot-start" initial conditions 
and those produced from dynamical models of the 
core- accretion process, based on the physics that 
they currently incorporate (J. Fortney 2010, pri- 
vate communication). 



2.3. Simulated Instrument 

Each target is observed in the J, H, and K 
bands with a hypothetical D — 8m telescope that 
uses a coronagraph. Companions may be detected 
exterior to the coronagraph inner- working- angle 
(IWA) which we assume to be 4 X/D. Corona- 
graphic occulting spots nominally have a radius of 
~ 3 X/D. In practice, however, the IWA is usually 
set by the type of speckle-suppression employed 
( Marois et al-lbood; iBiTler et al.lboioHCrepp et ail 



2010 



2011) 



We use an effective IWA of 4 X/D as 
a representative number. A case where the IWA 



is set to 3 X/D is considered in 

We assume that an "extreme" AO system 
generates a dark-hole search region centered on 
the s tar (jTrauger and Traubll2007t iBouchez et al. 
2009t ). The outer- working-angle (OWA), which 
is governed by the actuator-density of the high- 
order deformable mirror in the AO system, is set 
to 26 X/D. PSF subtraction is accounted for by 
extending the search area by several diffraction 
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widths from the OWA to 1.7" at the detector 
field-of-view (FOV) edge. It is assumed that the 
contrast remains constant in this region. The 
IWA and FOV edge set hard limits to the location 
where planets may be detected. 

Rather than explicitly modeling an optical sys- 
tem, we parameterize the contrast within the 
dark-hole such that sensitivity improves mono- 
tonically with angular separation from the star, 
falling by 1.5 mags fro m the IWA to OWA 
(jKataria and Simon l2010h . We define Co as the 
5cr contrast at the IWA in the J-band for a V = 5 
star, where a is the speckle noise standard devia- 
tion measured relative to the stellar peak intensity. 
When considering the H and K bands, the instru- 
ment IWA and OWA are scaled by (A/Ao) and 
the contrast within the search region is scaled by 
(A /A) 2 , where A = {1.25, 1.65, 2.20} /im is a ref- 
erence wavelength at the center of the band, such 
that longer wavelengths provid e wider and deeper 



is p 

dark-holes (|Malbet et ai1[l995l ) 



Contrast levels are further parameterized based 
on the target visual magnitude to account for the 
effects of AO system bandwidth, since fainter stars 
require longer wave-front sensor integration times 
which degrades sensitivity. The intensity of scat- 
tered starlight is scaled as the square of the rms 
wave -front error ( Malbet et ah[ 1995t Crepp et al 



2009) using calculations from lBaraned ([20081 ) that 
include the effects of time delay. The peak planet 
intensity in each band is likewise modified to ac- 
count for reductions in Strehl ratio as a function of 
target brightness. We assume that the AO system 
has an ability to optimize the wave-front sensor 
spatial sampling based on the target visual mag- 
nitude and assert that targets with V < 5 each 
achieve the same level of correction. In a similar 
fashion, we also consider the case of a near-infrared 
wave-front sensor (NIRWFS) that operates in the 
J-band, providing access to more low-mass stars 



( Rigaut et al. 19921 ). We account for the different 



number of photons per second arriving in the V 
and J bands for individual targets. 

Planets are detected when their brightness at a 
given separation exceeds the local 5<r contrast in 
either one of the three bandsH Planet apparent 



2 The cond03 and fort08 models are inconsistent with one 
another in the J and K bands, so we do no t perfo rm a 
bandpass search optimization study (c.f. Agol (2007)). 



magnitudes must also exceed the noise floor set 
by the sky background, which we assume to be 
m Bky = 24 in each band unless otherwise noted, 
corresponding roughly to the faintest detectable 
source in wl hr of integra tion time using an 8- 1 0m 
class telescope with AO (jUchimoto et al.ll2008l ). 

When simulating next-generation instruments 
that use "extreme" AO, we consider contrast lev- 
els: Co = 5 x 10~ 6 , 5 x 10~ 7 . Most current sys- 
tems employ low actuator-density deformable mir- 
rors and are limited by non-common-path errors 
between the AO wave-front sensor and science in- 
strument, limiting contrast levels to Co = 10~ 3 — 
10~ 5 . We also perform runs using Co = 5 x 10~ 5 
and a wider FOV to compare our calculations to 
results from current observing programs. 

3. Results 

In the following we present the output of our 
Monte Carlo simulations for several different high- 
contrast imaging surveys. Figures generally dis- 
play two plots: the number of planets detected 
in each stellar mass bin, < dn p /dM* >, and the 
number of planets detected per star in each stellar 
mass bin, < (drip/dN^/dM* >, where brackets 
denote the ensemble average over ss 10 3 realiza- 
tions. For each survey, we consider different val- 
ues for a and /? from Eqn. [TJ maximal range of 
semimajor axes, power- law index for planet masses 
around A-stars, and observations using a NIR- 
WFS. Adjacent curves within each plot nominally 
have all but one variable held constant to quan- 
tify the relative impact of various effects on the 
number and efficiency of detections. 

3.1. Volume-Limited Survey 

Volume-limited surveys may be useful when a 
large number of stars are considered, because they 
sample a wide range of relevant stellar parame- 
ters (mass, age, metallicity) , and the results are 
easier to interpret and less sensitive to uncertain- 
ties in the measured properties for specific targets 
compared to intentionally biased surveys. Simula- 
tions of a high-contrast imaging survey of 10,081 
A1-M2 V-IV stars within 50 pc (see JO for de- 
tails) are shown in Figs. 2, 3, using a contrast of 
Co = 5 x 10~ 7 with the cond03 and fort08 planet 
evolutionary models respectively. 
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3.1.1. Planet Detections 

We find that the number of directly detectable 
gas-giant planets in the solar neighborhood fol- 
lows a distribution that is multi-modal when plot- 
ted against stellar mass, with peaks located at 
M* < O.4M , M* ps l.OM , and M* ps 1.6M q . 
The middle peak is common to all cases consid- 
ered, while the last peak occurs only when taking 
into consideration correlations between star mass 
and planet mass. The first peak in the number of 
planet detections is seen in cases where the planet 
occurrence rate scales with star mass (Equ. 1), 
while its existence is only implicit in other situa- 
tions (see Fig. 2). 

These results are driven by the tradeoffs listed 
in SJTJ For late-type stars, the number of detec- 
tions is governed by the competing effects of an 
increasing stellar mass function, favorable star-to- 
planet brightness ratios, a diminished planet oc- 
currence rate, and poor AO correction. To provide 
a fiducial measure for the number of available tar- 



gets, we have over-plotted the iReid et al.l ((2002) 
PDMF, scaled down by a large multiplicative fac- 
tor to fit inside of Figs. 2, 3. The number of de- 
tections closely follows this shape for early M-type 
and late K-type stars. In this regime, targets are 
sufficiently faint such that nearly all young Jovian 
planets in wide orbits have an excellent chance of 
being imaged directly, provided they are located 
inside of the instrument FOV. The factor that lim- 
its the number of detections is that few nearby 
late-type stars are also youngH 

The number of planet detections begins to turn 
over (decrease) for M-stars when the planet oc- 
currence rate obeys Equ. 1, owing to the fact that 
low-mass stars have few Jovian planets. Further, 
insufficient compensation for atmospheric turbu- 
lence for the faintest stars reduces the Strehl ra- 
tio and can lead to prohibitively small companion 
peak intensities when observed against the NIR 
sky background. As a result of these effects, the 
number of detectable planets must approach zero 
at M, < 0.4M Q , even when using a NIRWFS. 
However, depending on the formation scenario, for 
instance whether or not a and /3 have non-zero 
values, a peak may not occur until much lower 
masses, possibly in the brown dwarf regime. 



3 See lLiu et al.l II2009T ) for a discussion of the "missing young 
M-dwarfs problem". 



Solar-type stars are sufficiently bright for "ex- 
treme" AO correction, but the number of available 
targets decreases sharply for M* > l.OM where 
the PDMF in the solar neighborhood experiences 
a break in the power law index. A local maxi- 
mum in the number of planet detections occurs 
at ps 1.0 Mq independent of simulation input. 
The number of planet detections falls monotoni- 
cally for M* > 1.0M Q , except when there exist 
correlations between star mass and planet mass. 
We find that changing the index of the planet 
mass power law from 7 a = — 1.4 to 7.4 = —1.0 
for M* > 1.5Mq generates another peak near 
M* ~ I.6M0 by amplifying the number of de- 
tections around massive stars. 

Increasing the range of semimajor axes from 
a = 35 AU to a — 120 AU enhances the num- 
ber of planet detections for most stellar masses. 
However, the most massive stars (M» ps 2.5M ) 
may experience the smallest relative gain since 
we have scaled the maximal extent of semima- 
jor axes by star mass (|}2]). The instrument FOV 
edge is 1.7" from the star and planets with a = 
120 AU and near face-on orbits around A-stars 
within 50 pc will subtend a wide angle on the sky 
and lie outside of the search region. This effect, 
though minor, will also impact next generation 
high-contrast imaging instruments wh en observing 
stars with close proximity to the Sun (jCrepp et al. 
2009). For example, the GPI has a FOV width of 
2.8" x 2.8" and will therefore not be able to detect 
and characterize Fomalhaut b in default operating 
mode, i.e., with the star occulted by the corona- 
graph. 

We find that a NIRWFS operating in the 
J-band can increase the number of planet de- 
tections around low-mass stars by tens of per- 
cent, providing access to more targets with red 
spectral energy distributions. Presumably, a 
laser-guide-star system can serv e a similar role 
(jWizinowich and Campbe 11 2009). Stars with blue 
spectral energy distributions experience only a 
modest loss. We note that such considerations 
are less important for space observations because 
the stellar wavefront may be sensed in the image 
plane b y the science camer a at the same wave- 
lengths (jGive'on et al.ll2007t ). 

The relative number of detections as a function 
of stellar mass is qualitatively similar between the 
top panels of Fig. 2 and Fig. 3, but the abso- 
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Fig. 2. — Monte Carlo simulation results of a 50 pc groun d-based volume- limited survey of 10,081 stars using 
the planet occurrence-rates from I Johnson et al.l (|2010al ). stellar and planetary distributions in Table [IJ an 
8m telescope with "extreme" AO and a NIR coronagraph that generate contrast levels of Co = 5 x 10 -7 , and 
the cond03 planet thermal evolutionary models. The overall number of detectable planets (top) is governed 
primarily by the stellar mass function (PDMF) for late spectral types, while observing efficiency, or the 
number of planet detections per star (bottom) , lies heavily in favor of early spectral types as a consequence 
of their intrinsic youth and visual brightness for AO correction. Including correlations between star mass and 
planet occurrence rate and other properties exacerbates the difference. Legends are the same for both plots. 
Main-sequence spectral-types are labeled on the top horizontal axis for reference. Each curve is sampled by 
21 data points (histogram bins) in stellar mass. 



lute numbers differ because the cond03 and fort08 
models make different estimates for the bright- 
ness of young planets. While the models generally 
agree on timescales of several Gyrs, the brightness 
difference can be an order of magnitude for mas- 



sive planets with t ago < 100 Myr and two orders o f 
magnitude at t agc < 10 Myr (jMarlev et alJl2007t l. 
The disparity between model predictions is evi- 
dent for bright stars, where contrast matters most, 
and with instruments that generate contrast levels 
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Fig. 3. — Same as Fig. 2 using the fort08 evolutionary models. Fewer planets are detected compared to the 
"hot-start" models but the relative number and rate of detections are qualitatively similar. 



of Co > 5 x 10~ 6 . In this regime, where only the 
youngest and most massive planets are detectable, 
the cooling tracks are most discrepant, requiring, 
for example, approxi mately 1 Gyr to conv erge for 
an m p = 8Mj planet (|Fortnev et alj|2008ft . Table 
[2] shows the total number of planet detections for 
each case considered in different stellar mass bins, 
using C = 5 x 1CT 6 and C = 5 x 10~ 7 and two 
different scmimajor axis distributions. 



3.1.2. Planet Detection Rates 

Importantly, the shape of the planet detec- 
tion curves do not follow the PDMF at high star 
masses, but instead rise above the (normalized) 
stellar mass function for all cases considered, in- 
dicating significantly higher detection efficiencies 
compared to low-mass stars. Dividing the number 
of detections in each stellar mass bin by the num- 
ber of stars observed in each stellar mass bin to 
calculate a detection rate, < {dn p /dN*)/dM* >, 
we hnd that the likelihood for imaging a planet fa- 
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Table 2: Number of planet detections (ensemble average) for a 50 pc ground-based volume- limited survey of 
10,081 stars over several stellar mass bins. Entries for cond03 and fort08 models are shown in the left and 
right of each column respectively. Two different semimajor axis distributions are considered. The gas-giant 
planet occurrence rate, k e — {0.17, 0.28}, is applied to semimajor axis distributions truncated at a = 35 AU 
and a = 120 AU respectively. Calculations allow for multi-planet systems. Rows having parameters that 
may be considered as the baseline case for an extrapolation of the RV planet population are bold-faced. 



vors massive stars by more than an order of mag- 
nitude over low-mass stars. Detection rate curves 
are shown on a logarithmic scale in the bottom 
panels of Fig. 2, 3. 

Massive stars have a high planet occurrence- 
rate and are comparatively young while on the 
main-sequence and early subgiant branch. Given 
their youth, they have bright planets and are also 
slightly more metal-rich than low-mass stars (on 



average), which leads to an additional increase 
(coupling) in the planet occurrence-rate through 
Eqn. [TJ Massive stars are also expected to have 
wide planet formation zones and are sufficiently 
bright at visible wavelengths to serve as their own 
guide star. As stars with progressively higher 
masses are observed, these effects combine in a 
non-linear way to increase the efficiency of direct 
imaging observations. As a result, A-stars to yield 
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a comparable number of planet detections as M- 
stars despite comprising a much smaller fraction 
of potential targets in the sky. 

While volume limited surveys may at first seem 
appealing because they involve a systematic search 
of the Sun's closest neighbors, they are inefficient 
unless targets are restricted to high stellar mass. 
For instance, using the cond03 models with Cq = 
5 x 1CT 7 , a = 120 AU, ~f A = -1.0, a = 1.0, and 
j3 = 1.2, the average number of A/* = 1.5M Q stars 
within 50 pc that need to be observed to discover 
one Jovian planet in an unbiased survey is «15, 
while the corresponding number at M* = 0.7M Q 
is «85 with a NIRWFS, and «107 without. In 
i j3.2l and ij3.5l we consider surveys where stars are 
pre-selected for youth to make a more interesting 
comparison and determine whether this trend with 
stellar mass is upheld. 

3.2. Age and Brightness Selected Survey 

Most current and planned surveys select tar- 
gets to maximize the chances for discovery. In 
this section, we consider another large-scale high- 
contrast imaging program in which stars from the 
50 pc volume-limited sample are chosen to have 
apparent magnitudes, V < 9 (J < 8 with a NIR- 
WFS), and ages, t agc < 1 Gyr. An age cut of 1 Gyr 
is convenient because it corresponds to a value in 
which discrimination between marginally younger 
or older stars may be done with reasonable consis- 
ten cy based on their l evel of chromospheric activ- 
ity (|Noves et al.| [l984f). Further, low -mass planets 
that are 1 Gyr old have likely cooled to a point 
beyond which direct detection becomes extremely 
challenging. Selecting bright and young stars from 
the volume- limited sample yields a total of 610 tar- 
gets when using a VISWFS and 760 targets when 
using a NIRWFS. Our results for simulated obser- 
vations of each star with Co = 5 x 10~ 7 using the 
same techniques and parameters described in 
and $3.1\ are shown in Figs. 4, 5. Table 3 dis- 
plays the number of planet detections for various 
contrast levels and two different semimajor axis 
distributions. 

3.2.1. Planet Detections 

We find that the number of planet detections 
rises steadily with increasing stellar mass from M- 
dwarfs to G-dwarfs and, like a volume- limited sur- 



vey, exhibits a peak near l.OM . However, unlike 
a volume-limited survey, few detections are made 
at low stellar masses on account of a dearth of 
available targets. This result is a consequence of 
the local age distribution within the solar neigh- 
borhood (Fig. 1). Only < 15% of nearby FGK 
stars are younger than 1 Gyr, and M-st ars are 
older yet (jMamaiek and Hillenbrandll2008l ). 

Of the propitious K and M-stars that are gen- 
erated in realizations of our age and brightness- 
selected samples - those that are young, nearby, 
and have a planet (e.g., analogous to e Eridani), 
only a fraction will have orbit orientations that 
are conducive to direct detection. While other 
spectral-types experience similar geometric losses, 
because many planets are seen in near edge-on 
orbits and remain hidden behind the corona- 
graphic mask (e.g., (3 Pic b), this final reduc- 
tion is enough to essentially preclude planet de- 
tections around M-stars unless a NIRWFS is em- 
ployed. To illustrate the difficulty in finding viable 
low-mass targets, we count those available from 
IZuckerman and Song! ( 2004 ) and IShkolnik et al 



(|2009l ) . The total number of single K and M-stars 
with V < 9 and s < 50 pc, including the /3 Pic- 
toris, TW Hydrae, Tucana-Horologium, r\ Cha, 
and AB Dor groups, is only 16. With planet oc- 
currence rates < 10% at wide separations, it is 
difficult to make a substantial number of detec- 
tions with such small samples. The GPI will use 
a wavefront sensoi@ that operates in the I-band 
and can therefore expect a number of planet de- 
tections that lies in between the (brown dash and 
red dash-dot) curves shown in Figs. 4, 5. 

A scenario in which correlations between star 
mass and planet occurrence rate are weak can 
improve the number of detections at low star 
masses. Raising the overall planet occurrence rate 
can also help, although the size of the error-bar 
on k from Eqn. [T] is already small. Previous 
high-contrast imaging observations have placed 
firm upper-limits on the overall frequency of mas- 
sive (m„ > 4Mj) planet s in wide orbits at 20% 
(jNielsen and Close! l2010h . Otherwise, loosening 



4 Most "extreme" AO systems will employ two wavefront sen- 
sors: one that is sensitive to dynamic wavefront errors and 
runs at a frequency /ail, 000 Hz, and another that is more 
sensitive to static er rors, running at a frequency of /~1 Hz 
(Wallace et al. 2010). Here we arc referring to the dynamic 
wavefront sensor. 



12 




Fig. 4. — Monte Carlo simulation results of a 50 pc survey of 610 st ars (760 with NIRWF S) with age, 
i ago < 1 Gyr and V < 9 (J < 8) using the planet occurrence rates from Johnson et al.l ( 2010al ). stellar and 
planetary distributions from Table 1, an 8m telescope with "extreme " AO and a near- infrared coronagraph 



that generate contrast levels of Co = 5 x 10 , and the iBaraffe et al.l (J2003J) "hot-start" planet evolutionary 
models. A-stars yield a comparable number of detections to FGKM stars and have elevated detection rates 
even when stars are pre-selected for youth and brightness. Legends are the same for both plots. Main- 
sequence spectral-types are labeled on the top horizontal axis for reference. Each curve is sampled by 21 
data points in stellar mass. Color-scheme is the same as the volume-limited case. 



age and brightness selection cuts can enhance the 
number of planets imaged around low-mass stars 
at the expense of detection efficiency. Results from 
the volume-limited case (Fig. 2, 3) show the num- 
ber and rate of detections in the limit as both cri- 
teria are completely removed. 

The distribution of planet detections around 



massive stars is similar to that of a volume-limited 
survey, indicating that a significant fraction of 
stars with M* > 1.8M have made the cut be- 
cause they are intrinsically young and bright, ow- 
ing to the applied selection bias. Fig. 6 shows 
the relationship between the median and maxi- 
mum age of stars used for our simulations. Given 
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Fig. 5. — Same as Fig. 4 using the Fortney et al. 2008 "core-accretion" models. 



the difficulty in age-dating stars in the field, stellar 
mass serves as an excellent proxy and can aid in 
the selection of promising high-contrast imaging 
targets for a number of reasons. 

3.2.2. Planet Detection Rates 

Dividing the number of planet detections in 
each stellar mass bin by the number of targets ob- 
served in each stellar mass bin, we find that the 
probability for imaging a planet is enhanced for 
all spectral-types compared to the volume-limited 
case (as expected), and that massive stars are 
again preferred targets (Fig. 4, 5). Both mod- 
els indicate however that this result is sensitive to 



the value of a, the power-law index that controls 
the planet occurrence rate as a function of stellar 
mass. When a = 0, the detection rate becomes 
relatively flat in the range, 1.0 < M»/Mq < 2.6, 
and exhibits a small negative slope (i.e., efficien- 
cies favoring lower mass stars) for M» < I.OMq. 
Although the slope changes sign in a regime where 
few planets are imaged, a NIRWFS can provide ac- 
cess to a sufficient number of stars to notice such 
a trend in practice when Co ~ 5 x 10 ~ 7 . Corre- 
lations between metallicity and planet occurrence 
rate can further improve number statistics at low 
star masses. Dropping our assumption that semi- 
major axis extent scales with stellar mass yields 
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1.0 


1.2 


5 x 10~ 7 


-1.0 


120 


nir 


18 


19 


30 


26 


24 


15 


12 


5.7 


dn p /da oc a u , k e = {0.17,0.27} 


input parameters 




stellar 


mass 


range 


(M Q ) 




a 


P 


Co 


lA 


a 


WFS 


0.40-0.95 


0.95-1.50 


1.50-2.05 


2.05-2.60 








5 x 10" b 


-1.4 


35 


vis 


0.7 


0.1 


1.3 


0.0 


0.5 


0.0 


0.2 


0.0 





1.2 


5 x 10~ 5 


-1.4 


35 


vis 


1.1 


0.2 


2.1 


0.0 


0.8 


0.0 


0.4 


0.0 


1.0 


1.2 


5 x 10~ 5 


-1.4 


35 


vis 


0.8 


0.2 


2.4 


0.0 


1.4 


0.0 


0.9 


0.0 


1.0 


1.2 


5 x 10" 5 


-1.4 


120 


vis 


2.7 


0.5 


9.3 


0.3 


5.1 


0.0 


2.7 


0.0 


1.0 


1.2 


5 X 10~ 5 


-1.0 


120 


vis 


2.7 


0.5 


9.3 


0.3 


8.7 


0.0 


4.6 


0.0 


1.0 


1.2 


5 x 10" 5 


-1.0 


120 


nir 


6.9 


2.1 


8.1 


0.2 


7.1 


0.0 


4.1 


0.0 








5 x 10" B 


-1.4 


35 


vis 


1.7 


0.7 


3.5 


0.3 


1.2 


0.0 


0.5 


0.0 





1.2 


5 x 10~ 6 


-1.4 


35 


vis 


2.7 


1.4 


6.0 


0.6 


2.3 


0.1 


1.0 


0.0 


1.0 


1.2 


5 x 10~ 6 


-1.4 


35 


vis 


2.2 


0.9 


7.0 


0.7 


4.0 


0.1 


2.3 


0.0 


1.0 


1.2 


5 x 10~ 6 


-1.4 


120 


vis 


3.7 


1.7 


16 


3.2 


8.0 


0.7 


4.2 


0.2 


1.0 


1.2 


5 X 10~ 6 


-1.0 


120 


vis 


3.7 


1.6 


16 


3.1 


13 


0.9 


6.5 


0.2 


1.0 


1.2 


5 x 10~ 6 


-1.0 


120 


nir 


11 


6.3 


15 


2.0 


10 


0.2 


5.6 


0.0 








5 x 10~ 7 


-1.4 


35 


vis 


3.7 


3.9 


10 


8.8 


4.1 


2.7 


1.6 


0.9 





1.2 


5 x 10~ 7 


-1.4 


35 


vis 


6.4 


6.8 


18 


15 


7.0 


4.8 


2.9 


1.6 


1.0 


1.2 


5 x 10~ 7 


-1.4 


35 


vis 


5.1 


5.4 


21 


18 


12 


8.3 


6.6 


3.6 


1.0 


1.2 


5 x 10~ 7 


-1.4 


120 


vis 


8.4 


8.7 


43 


41 


22 


18 


11 


7.5 


1.0 


1.2 


5 X 10~ 7 


-1.0 


120 


vis 


8.4 


8.7 


43 


41 


33 


27 


16 


11 


1.0 


1.2 


5 x 10~ 7 


-1.0 


120 


nir 


23 


25 


40 


36 


28 


19 


14 


7.4 


Number of Stars Observed with VISWFS 


87 


331 


139 


53 


Number of Stars Observed with NIRWFS 


228 


340 


139 


53 



Table 3: Number of planet detections for a 50 pc survey that targets 610 stars (760 with NIRWFS) with 
i ago < 1 Gyr, and V < 9 ( J < 8). Substantial gains are attained when contrast levels approach Co ~ 5 x 10 . 
Entries for cond03 and fort08 models are shown in the left and right of each column respectively. Two different 
semimajor axis distributions are considered. The gas-giant planet occurrence rate, k e — {0.17,0.28}, is 
applied to semimajor axis distributions truncated at a = 35 AU and a = 120 AU respectively. Calculations 
allow for multi-planet systems. Simulations with Co = 5 x 10~ 5 use a 200 X/D wide FOV. Rows having 
parameters considered as the baseline case for extrapolation of the RV planet population are bold-faced. 
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Monte Carlo simulations to ensure consistency 
with (pre-)main-sequence or subgiant branch evo- 
lution. Intrinsically young, massive stars have 
bright planets compared to low-mass stars, in- 
creasing the likelihood for a direct detection. In 
the absence of information regarding membership 
to a moving group or association, stellar mass may 
be used as a trivial aid in the selection of viable 
field stars for high-contrast imaging targets. 



slightly smaller (more negative) slopes (see §3.5|) . 

Given these results, we conclude that it will 
be important to carefully examine the detection 
rate of future high-contrast imaging programs as 
a function of stellar mass by making plots sim- 
ilar to Fig. 4, 5. Such diagnostics require a 
large sample of targets but can potentially help 
to discriminate between planet formation scenar- 
ios. When targets are chosen carefully, so as to 
minimize the effects that stellar age and bright- 
ness have on detection efficiencies, the dependence 
of the planet occurrence rate upon star mass (Equ. 
1) is strong enough to help assess whether the pro- 
cesses responsible for forming the RV planet pop- 
ulation (a < 6 ATI) may also operate at larger 
separations. 

Planet evolutionary models are sensitive to age, 
and stellar ages are not oriously difficult to de- 
termine with accuracy (ISoderbloml l2010h . We 
find that our results are robust to perturba- 
tions to the distribution in Fig. 1, provided that 
there remains some semblance of a roughly con- 



stant star-formation rate over t he past 4-6 Gyr 
(jMamaiek and Hirlenbrandll2008h . 



3.3. Current Technology 

Motivated by the direc t detection of planets 



orbiting f3 Pictori s (A6V) ([ Laarang e et al 



planets 

mmk, 



HR 8799 (A5V) dMarois et al. 2008, 2010), and 
FomalhaulO fA3V) (|Kalas et all 120081 ). we have 
performed additional calculations using contrast 
levels commensurate with those achieved by in- 
struments using standard AO, to compare our 
results with recent observations. Repeating the 
age-and-brightness selected survey where targets 
within 50 pc are chosen based on their youth, 
i agc < 1 Gyr, and brightness for AO correction, 
V < 9 (J < 8), Fig. 7 shows our results for 
the number and rate of planet detections using 
the cond03 models when setting contrast levels to 
Co = 5 x 10 -5 and expanding the search region 
to ±100 \/D (11// wide) from the star to bet- 
ter reflect the typical FOV of current instruments 
( Marois et al 1 120081: iBiller et al.ll2010h . The num- 
ber of planet detections using both planet evolu- 
tionary models and various input parameters is 
also shown in Table 3. 

We find that degrading contrast levels yields re- 
sults that resemble Figs. 4, 5, only the number and 
efficiency of detections is essentially scaled down 
by a factor of several for all curves. Careful exami- 
nation between the two cases shows that low-mass 
stars have gained some relative ground when Co 
is adjusted from 5 x 10~ 7 to 5 x 10~ 5 , on account 
of less demanding contrast requirements compared 
to massive stars, but nevertheless the results are 
qualitatively similar. Low-mass stars generate few 
detections, though are relatively efficient targets 
when a — 0; the detection rate curve is fiat be- 
yond M* = 0.9M© when a = 0; and high-mass 
stars yield a comparable number of detections as 
low-mass stars but with high detection rates when 
a = 1, particularly when the power-law index gov- 
erning planet masses is modified from 7^4 ~ —1.4 
to 7a = -1.0. 

In other words, setting Co = 5 x 10~ 5 we 
find that naive extrapolation of the known planet 
population, i.e., those discovered by the RV 
technique, to separations wider than several 



5 Fomalhaut b has mysterious colors and has yet to be de- 
tected from the ground in the NIR. 
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Fig. 7. — Planet detection numbers and rates as a function of stellar mass for a 50 pc age-and-brightness 
selected ground-based survey using contrast levels of Co = 5 x 10 -5 , comparable to present-day instruments. 
Low-mass stars yield few detections and high-mass stars offer the highest efficiency when planets at wide 
separations resemble the RV sample. For reference, 95% confidence levels from a number of high-contrast 
surveys yielding non-detections (V < 9, t ago < 1 Gyr) are over-plotted in the bottom panel. Constraints 
for M-stars represent a lower-limit since we have included V < 11 targets in this bin. Our calculations are 
thus consistent with high-contrast imaging surveys conducted to date in terms of both relative and absolute 
numbers. Simulation results may be scaled up or down based on the actual occurrence rate and truncation 
radius of planets in wide orbits (e.g., consider the difference between the solid blue and solid green curves). 



AU is consistent with A-stars dominating the 
first planet imaging discoveries (see also §3.5[) . 
This interpretation of our simulations is rein- 
forced by the null-results of other high-contrast 
progr ams that have targete d nearb y solar-type 
stars: iMasciadri etail (1 20051) / M05) iBiller et al ' 



(200 71) (B07). lLafreniere et all (120071) (L07), and 



iMetchev and Hillenbrand! ~ 20091 ) (deep sample 
only) (MH09) selected targets based primarily on 



proximity and age and were marginally sensitive 
to wide-separation planets, with Co « 10~ 4 -10 -5 . 
Combined, these surveys observed: 50 K-stars, 61 
G-stars, 16 F-stars, but only 1 A-star within 50 
pc {V < 9, i agc < 1 Gyr)H 

The combination of null results from these early 
high-contrast imaging surveys along with the de- 



3 Consult iNielsen and Closd j2010l) for more details. 
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Fig. 8. — Cumulative number of planets detected as a function of stellar mass for a 50 pc age-and-brightness 
selected survey using the cond03 models (top) and fort08 models (bottom). Several characteristic contrast 
levels are considered. The number of detections using current technology well-matches the results from obser- 
vational surveys to date, predicting between 0.6-16.6 planet imaging discoveries for Co = 5 x 10~ 5 . Surveys 
resembling the NICI campaign can at most double the number of planet imaging detections. Substantial 
gains occur when contrast levels improve beyond Co = 5 x 10~ 6 and low-mass planets («0.5Afj) are imaged. 



tection of massive planets with wide orbits around 
subgiant stars (former or "retired" A-stars) us- 
ing the RV technique has prompted other groups 
to focus mor e -so o n massive stars. Recently, 
Leconte et al.l ( 2010l ) (L10) and Ehrenreich et al 



(|2010h (E10) have published non-detections for 
substellar companions orbiting: 3 G-stars, 14 F- 
stars, 22 A-stars, and 4 B-stars within 60 pc 
(V < 9, i ago < 1 Gyr), achieving similar contrast 
levels. As shown in Fig. 7, the average number 
of young, nearby stars required to detect a large- 
separation planet, given our assumptions and no 



other previous knowledge on the presence of com- 
panions or a debris dislo, is higher than that re- 
ported for a < 120 AU (95% confidence level). 

Using parameter values a = 1, /3 = 1.2, a = 
120 AU, 7a = -1.0, and dn p /da oc a~ 61 with 
a VISWFS, the fort08 and cond03 models yield 
0.6 ± 0.2 and 16.6 ± 4.4 total planet detections 
(ensemble average), taking into account the uncer- 
tainty in parameters from Equ. 1 (fc e = 0.27±0.04, 



7 HR 8799, Pictoris, and Fomalhaut are clearly special 
cases in this regard, each having prominent debris disks. 
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a = 1.0 ± 0.3, /3 = 1.2 ± 0.2). The number 
of directly imaged planets orbiting young stars 
within 50 pc currently stands at 5 (HR 8799 
bcde, (3 Pic b). Our simulations using the "core- 
accretion" and "hot-start" thermal evolutionary 
models, which represent extrema in brightness 
predictions for young planets (see §2.2|) . are thus in 
excellent agreement with observational results to 
date in terms of both relative and absolute detec- 
tion rates, straddling the number of known plan- 
ets with expectations for several more discoveries 
in the near-term. 

Fig. 8 shows the cumulative number of planet 
detections as a function of stellar mass for various 
contrast levels using the cond03 models and same 
input values listed above. Both a VISWFS and 
NIRWFS are considered for next-generation in- 
struments. The Near-Infrared Coronagraphic Im- 
ager (NICI) campaign occupies an interesting pa- 
rameter space and is conducting a thorough search 
of ~ 300 young stars using the Gemini S. AO sys- 
tem, a coronagraph, and speckle suppression tech- 
niques that gene r ate unprecedented contrast levels 



(ILiu et all 120101: iBiller et al.ll2010t IWahhaj et al 



20111 ). Despite these advances in hardware and 
data-processing, routinely achieving sensitivity to 
> lMj companions, our simulations indicate that 
NICI and other surveys of similar scope will at 
most double the number of planet images. There 
exist solar-type stars within 50 pc with imageable 
planets at Co ~ 5 x 10~ 6 , but the efficiency of ob- 
servations in this regime is low, requiring a dedi- 
cated effort involving observations of all available 
nearby young FGK stars, particularly if planet 
semimajor axes extend only to several tens of AU. 

We predict that substantial gains will be ob- 
tained only when contrast levels improve by two 
orders of magnitude compared to observations us- 
ing standard AO. Instruments that employ "ex- 
treme" AO systems will increase the number of 
planets imaged directly around V < 9 (J < 8) 
stars within 50 pc, yielding between 12-81 discov- 
eries total. Observations of fainter stars can fur- 
ther boost the number of detections at the expense 
of many more nights at the observatory (£ !3.1j) . 

3.4. More Distant Stars 

Several dozen extrasolar planets orbiting stars 
within 50 pc of the Sun will be imaged and char- 
acterized by high-contrast instruments over the 



next few years. While preliminary trends may 
begin to emerge f rom these datasets (e.g., see 
iMarcv et al. (1999) for a review of the RV pop- 
ulation when it comprised 17 exoplanets), it is im- 
portant to collect the largest possible sample for 
statistical analyses to inform our understanding of 
gas-giant planet formation and evolution. Indeed, 
planet properties and correlations with host star 
properties are likely to involve distributions that 
are multi-dimensional, multi-modal, and possibly 
time-dependent at the youngest ages. 

We have performed calculations for young and 
bright stars beyond 50 pc to gain a handle on the 
expected yield of magnitude-limited high-contrast 
surveys. Results for the planet detection rate 
setting Co = 5 x 10~ 7 are shown in Fig. 9 
for stars separated into three different mass bins, 
1.0 < M»/Mq < 1.5, 1.5 < M*/M Q < 2.0, and 
2.0 < M*/Mq < 2.5, along with the cumulative 
number of detections as the survey maximum dis- 
tance is increased outwards in 10 pc increments. 
Targets are again chosen such that V < 9 ( J < 8 
with NIRWFS) and £ agc < 1 Gyr using the distri- 
butions from Table 1. Simulation input is selected 
to represent a "pessimistic" case in order to place 
a lower-limit on the number of planet detections 
gained. We have assumed that the semimajor axis 
distribution of the RV planet population extends 
only to a = 35 AU, and that the sky background 
is bright, limiting planet detections to apparent 
magnitudes 77ij HK <C 22. Further, we have removed 
the assumption that more massive planets orbit 
stars with mass M* > 1.5A/Q. 

Despite assuming non-ideal input parameters, 
we find that additional wide-separation compan- 
ions may be imaged around 1.0 < M*/M Q < 2.5 
stars well beyond 50 pc, the canonical distance 
that most planned surveys are restricted. Consis- 
tent with previous sections, the most massive stars 
(M* > 2.0M Q ) offer the highest detection rates. 
Average efficiencies for each group fall slowly with 
increasing distance and asymptotically approach a 
minimum value as the maximum number of avail- 
able stars in each mass bin is reached from the 
applied V < 9 brightness selection cut. The num- 
ber of detections does not plateau until distances 
exceeding 100 pc. 

Stars with mass, 2.0 < M*/M < 2.5, offer de- 
tection efficiencies of at least one planet detection 
per 17 targets out to a distance of 87 pc (aver- 



19 



O.lOr 



0.08 



n 

+J 
w 



£~ 0.06: 
Oft 

0.041 

Q V 

+J 

2 0.02 



V < 9 
age < 1 Gyr 



'*/*,. 1 detection per 17 stars ^^^^^^ 15 7 

**»,;, 



■1.0 < M,/ M -m , < 1.5 
1.5 < M */ M„„„ < 2 . 
■2.0 < M */ M, m < 2 .5 



0.0C 



cond03 
C =5 

dn /da 0= a 
a J 35 AO 

If -1 - 4 

I 

60 



is £^ 

"A... 23 4 ^ 

•.... 2 °- 7 27.6 -" 

"•At,... 1 28.8 "•'•Hi ,_ - 

w T; s-. 7 . ...... 



70 80 90 100 110 

survey maximal distance (pc) 



120 130 



Fig. 9. — Monte Carlo simulations of an age-and-brightness limited high-contrast imaging survey targeting 
stars with incrementally larger maximum distances. Such observations can provide an additional factor of 
at least 2 — 3 in the total number of planet detections for M» > 1.0M Q stars compared to canonical surveys 
restricted to 50 pc. Massive stars yield the best observing efficiency, offering more than one detection per 17 
targets for stars out to a distance of w 87 pc averaged over the entire volume. Detection rates approach a 
minimum as the survey becomes magnitude-limited due to the applied V < 9 target selection cut. The total 
cumulative number of detections is shown above each data point. For example, the number of detections 
between 60 — 70 pc for 1.5 < M*/Af Q < 2.0. stars is 15.5 — 11.1 = 4.4 planets. Simulations assume a sky 
background limit of 22 mag, semimajor axis distribution dn p /da oc a , with a < 35 AU extent, and 
planet mass power- law index 7^4 = —1.4. 



aged over the entire volume). For reference, ap- 
proximately 10 stars per good-weather-night may 
be observed with an 8-10m telescope using the 
ADI technique. The number of planets detected 
in each stellar mass bin increases from 11.3, 7.6, 
and 4.4 at 50 pc to 25.3, 19.2, and 12.5 at 80 pc 
respectively, suggesting that a factor of 2-3 gain 
in the wide-separation planet population is pos- 
sible by including more distant targets. We find 
that improving the IWA from 4 X/D to 3 X/D 
can improve the number of detections by another 
~11% in each mass bin over a similar space vol- 
ume. Although proper motions in this regime are 
small, requiring long time baselines for follow-up 
astrometrjffl, augmenting target lists by including 



Candidates have been confirmed as common proper- 
motion companions ar ound stars as distant as 160 pc 
l lLafreniere et al.ll20lfl) . Further, it may be possible to es- 



stars with distances between 50 pc and the nearest 
star-forming regions at «140 pc will maximize the 
overall science output of future instruments. 

3.5. A Nearby Association 

Members of nearby moving groups and associ- 
ations fall under the category of target selection 
by age because the constituent members are likely 
to have formed around the same time. When a 
cluster of young stars are also situated at roughly 
the same distance from the Sun, they represent 
a special case in high-contrast imaging and con- 
venient control for comparing the number and 
rate of planet detections to the volume-limited 
and age-and-brightness-selected surveys presented 
previously. In this section, we further elucidate 



tablish association indirectly using spectroscopic character- 
ization. 



20 



our results by simulating observations of a col- 
lection of nearby, young stars, each with an age, 
t agc = 30 ± 10 Myr, at a distance, s = 48 ± 7 
pc. Such a consortium may be considered as rep- 
resentative of the T ucana-Horologium association 
jrres et aI1l2008l) . 

We consider three diffe rent cases using a = 35 



AU, C = 5 x 10~ 7 , and a lMiller and Scald (Il979h 
IMF: 

(I) a =1,(3 = 1.2, a max = (M*/2.5M ) 4 / 9 a 

(II) a = 0,13 = 0, a max = (A/„/2.5A/ Q ) 4 / 9 a 

(III) a = 0, /3 = 0, a max = (M*/2.5M )° 5 

The first two cases were treated in the 50 pc 
volume-limited and age-and-brightness-selected 
surveys. In the third case, we relax our assump- 
tion that the planet formation zone grows with 
stellar mass by setting the maximal extent of 
semimajor axes a max = 35 AU for all targets. 
For each scenario, the above age and distance 
ranges are used as limit values for uniform sam- 
pling. Tuc ana-Horologium currently h as 49 known 
members ( Zuckerman and SonsJ [20041 ) . We simu- 
late observations of all members (V < 12.2) and 
construct 10 ensembles, a larger number than in 
previous sections, to reduce the statistical noise 
associated with the small number of targets. The 
met allicity of each s tar is set to [Fe/H] = —0.25 
dex ()Makarovll20"07b[) . Only one planet evolution- 
ary model is considered, because different models 
yield scaled versions of one another when the age is 
fixed; i.e., the number of planets detected for dif- 
ferent models is related by a multiplicative factor. 
Fig. 10 shows our results assuming "hot-start" 
initial conditions. 

We find that low-mass stars yield the most 
planet detections and that there is a systematic 
decline with increasing stellar mass. This simple 
dependence is a result of the features held in com- 
mon with all stars in the association. With the 
same approximate age and distance, the number 
of planet detections is governed primarily by the 
stellar IMF. Setting a = and f3 = changes the 
slope of the planet detection curve further in fa- 
vor of low-mass stars. The steepest slope occurs 
when a = 0, (3 = 0, and a max = (M*/2.5M Q )° a. 
In this case special case, the ratio of the number 
of planet detections for high-mass stars to that of 



low-mass stars is close to the ratio of stars from the 
IMF - for instance, n p (2.OM©)/n p (l.OM ) = 0.24 
compared to (2.0M Q /1.0M©)- 2 = 0.25 - though 
slightly smaller because massive stars command a 
deeper contrast ratio. This is effectively a verifi- 
cation of our code. Indeed, controlling for stellar 
brightness and "turning off' the effects from the 
AO system produces a slope in the number of de- 
tections that matches exactly that dictated by the 
IMF. 

Dividing the number of planet detections by the 
number of targets observed in each stellar mass 
bin, we find that massive stars are preferred tar- 
gets when occurrence-rate and semimajor axis cor- 
relations are applied. However, setting a = 
(changing (3 only creates a constant offset) and 
eliminating the correlation between a max and M* 
results in an average slope of —0.07 dex / M over 
the stellar mass range considered. This simplest 
case is the only one, including simulations from 
all previous sections, where M-stars both dominate 
the number of detections and have the highest ob- 
servational efficiency. 

The vast majority of Tucana-Horologium mem- 
bers have been searched previously with AO and 
a coronagraph. There is one companion with a 
mass at the planet/BD boundary that was de- 
tected dir ectly orbiting the K 2V star AB Pic 
(V = 9.2) (|Chauvin et al.ll2005l ). Performing sim- 
ulations that mimic present-day observations, by 
degrading contrast levels to Cq = 5 x 10~ 5 , in- 
creasing the FOV (see £13.31 for details) , and also 
setting a = 300 AU, we find that the three cases 
yield a cumulative number of 0.3 ± 0.1, 1.2 ± 0.1, 
and 1.5±0.1 imaged planets over the MK spectral- 
type range respectively. 

The Tucana-Horologium association is cur- 
rently incomplete at low stellar masses, because 
the faintest members have yet to be identified. 
One might therefore cautiously speculate that the 
detection of a single planet candidate thus far fa- 
vors case III. However, our simulations yield the 
same answer for cases II, III to within round- 
ing errors. Also, a different formation mechanism 
could produce a different overall planet occurrence 
rate and we have set k e = 0.17 for each case. Fur- 
ther, planet evolutionary models are too uncertain 
at this age to justify using calculations involving 
absolute numbers alone to discriminate between 
possible formation scenarios. 
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Fig. 10. — Monte Carlo simulations of a nearby young stellar association synthesized to resemble Tucana- 
Horologium, using Co — 5 x 10~ 7 , a = 35 AU, and the cond03 planet evolutionary models. Gas-giant planet 
occurrence-rates as a function of stellar mass, a, and metallicity, j3, are changed, along with the scaling of 
the size of the planet formation zone with stellar mass. Low-mass stars dominate the number of detections 
as a result of the common age and distance shared by constituent members. The detection rate rises towards 
higher stellar masses when correlations between star mass and planet properties are considered. The slope 
of detection rate curve is negative when assumptions about planet occurrence rate are dropped and we relax 
our assertion that the maximal extent of planet orbits scales with stellar mass. In this case, low-mass stars 
can yield the most detections and highest detection efficiencies. The value of detection rate slopes may be 
used as a diagnostic to help discriminate between planet formation scenarios. 



These results are consistent with the recent 
detection of large-separation, planetary-mass ob- 
jects found orbiting late-type members of mov- 
ing groups, associations, and star-forming regions. 
Table 4 shows a compilation of directly imaged 
planets (m p < 15Mj) to date according to the 
exoplanet encyclopedia (http://exoplanet.eu/). 
There currently exists a dichotomy whereby the 



host stars fall cleanly into two categories of nearby 
massive stars on one hand, and low-mass stars in 
more distant associations/clusters on the other 
hand. 

The preference for planet detections around 
low-mass stars may also be explained by an obser- 
vational bias resulting from small number statis- 
tics. For a given tightly bound cluster, there are 
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Star 


SpTy 


A/*/M Q 


m p /Mj 


V 


Assoc. 


Age (Myr) 


a (AU) 


s (pc) 


p/d 


Fomalhaut 


A3 


2.1 


<3 


1.2 


Castor 


200 


119 


7.7 


0.5 


P Pic 


AG 


1.8 


8± 5 2 


3.9 


P Pic 


10-12 


8-15 


19.3 


0.7 


HR 8799 


A5 


1.6 


5-10 


6.0 


Columba? 


20-150 


15-68 


39.4 


1.8 


AB Pic 


K2 


0.8 


13.5 ±0.5 


9.2 


Tuc-Hor 


30 


260 


47 


0.2 


GQ Lup 


K7 


0.7 


21.5 ±20.5 


11.4 


Lupus 


1 


103 


140 


0.2 


1RXS J1609 


K7 


0.7 


«8 


7=11.0 


USco 


2-5 


330 


145 


0.1 


GSC 06214 


Ml 


0.6 


13.5 ±2.5 


7=11.1 


USco 


2-5 


320 


145 


0.1 



Table 4: Directly imaged planetary-mass (m p < 15Mj) companions with stellar primaries as of February 25, 
2011 (http:/ /exoplanet.eu/). There currently exists an apparent dichotomy between the detection of planets 
around A-stars in the solar neighborhood and the detection of ultra- wide separation planets around late-type 
stars in clusters. This division may be roughly parameterized according to the ratio of the cluster size to the 
average distance of its members from the Sun, p/d. Our calculations show that massive stars are preferred 
targets when searching within the solar neighborhood (s < 50 pc), provided strong correlations between 
star mass and planet properties are present. Young stellar clusters represent a special case in high-contrast 
imaging and also suffer from small number statistics, yielding the most planet detections around MK stars, 
irrespective of such correlations. If ultra- wide separation planets form by a separate process than those 
orbiting close to their star, the planet detection rate for clusters will have a negative or approximately flat 
slope when plotted against stellar mass. This diagnostic may be used to he lp differentiate b e tween formation 
scenari o s. Numbers used for ages, d istances, and spatial scatter are from iBeichman et al. ( 201dh . IMakarov 
( 2007a ). Barrado v Navascued (|l998h . the Extrasolar Planets Encyclopedia, and NStED. 



generally too few members to have a reasonable 
chance of imaging a planet around a high-mass 
star, despite correlations that may exist between 
star mass and planet properties. For example, 
Tucana-Horologium is one of the largest nearby as- 
sociations and contains only 4 A-stars. This bias 
is exacerbated by the fact that: (1) half of cluster 
members are binaries, and high-contrast programs 
nominally avoid binaries; (2) most stellar clus- 
ters are located at 140 pc, effectively reducing the 
overall planet occurrence rate to its value at large 
orbital separations; (3) current estimates of the 
frequency of ultra-wide separation planetary-mass 
companions are sal-4% (Metchev and Hillenbrand 



comp anions are ~i -47o (|J> 
20091 llreland et al.ll201ll 



We find that the two host star groups, late-type 
and early-type, may be roughly separated accord- 
ing to p/d, the ratio of the characteristic cluster 
size (scatter) to the mean distance of cluster mem- 
bers from the Sun. When p/d is large and nearby 
targets comprised of different moving groups, asso- 
ciations, recently dispersed clusters, and stars that 
have formed in relative isolation - each having dif- 
ferent ages - are spatially mixed (often labeled as 
field stars beyond a given age or when exhibiting 
space motion that is unreconcilable with known 



kinematic groups), the resulting number and rate 
of detections follow the trends in Fig. 4, 5, and 7. 
When p/d is small and a group of tightly clustered 
stars that unambiguously share the same age and 
distance are observed, the resulting number and 
rate of detections follow the trends in Fig. 10. 

A change in the sign of the detection rate 
slope is a specific prediction that can be used 
to help differentiate between formation scenarios 
(Fig. 10). If a negative detection rate slope is 
found for wide separation planets orbiting MK 
stars in tight clusters, then these planets could not 
have formed in the same way as HR 8799 bcde 
or j3 Pic b, because we find that strong corre- 
lations between star mass and planet properties 
are required to account for their existence ( §3.31) . 
While measurement of a positive detection rate 
slope is an inconclusive result by itself, measure- 
ment of the slope value at the 10% level can offer 
a more definitive discriminatory test. Diagnostics 
based on stellar mass can complement other stud- 
ies, such as those involving host star metallicities, 
semimajor axis distributions, and inferences made 
from absolute detection rates. Several teams are 
now con ducting dedicated su r veys of star-forming 
regions (|Todorov et alj|20ld : llreland et al]|201ll 
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Lafreniere et al.ll201l[ ). Their results can help de- 
termine whether multiple formation channels are 
at work. 



at the tens of percent level on the stellar mass- 
dependent binary fraction and whether the indi- 
vidual components of wide binaries are targeted. 



4. Caveats 

In this section, we list important caveats and 
assumptions that have gone into the paper and 
describe their qualitative influence on the results: 

(1) While this study is concerned primarily with 
calculations of the relative behavior in the number 
and efficiency of planet detections as a function of 
stellar mass, we have also made claims about the 
absolute number of detections. Figures 2-5, 7-10, 
and Tables 2, 3 may under-estimate the number 
of possible planet detections for several reasons: 
contrast levels below Co = 5 x 10 ~ 7 may be pos- 
sible for bright stars; the frequency of multiple 
planet systems may be greater than 28%; systems 
with more than 4 Jovian planets may exist; our 
AO system model does not account for reduced 
differential chromatic aberrat ions when using a 
NIRWFS (|Guvon et al.l 120061 ): instruments will 



also use t he polarization signatu re of planets for 



detection ( Roelfsema et al. 2010l ); and, more than 
one planet formation mechanism seems likely. At 
the same time, we may have over-estimated the 
number of planet detections: a dn p /da oc a~ 61 
semimajor axis distribution may be too optimistic 
for wide-separation companions; it is doubtful that 
surveys will have sufficient time allocations to ob- 
serve several hundred stars in each of the JHK 
filters; and, depending on the spectral-type, stars 
may require 2-3 hours of integration time each to 
reach contrast levels of order Co = 10~ 7 , there- 
fore also limiting the number of targets observed. 
Further, the effects of orbit migration and planet- 
planet scattering can also alter the number of 
planet detections. 

(2) We have applied a uniform multiplicity frac- 
tion to stars of all masses. Low-mass stars appear 
to have a binary fraction that is 10-20% lower 
than solar-type stars as a consequence of natal 
gravitational interactions with cluster members 
((Fischer and Marcvl fl992l) . The binary fraction 
of high-mass stars is plausibly h igher than 0.33, 



(3) We have calculated the number of planet de- 
tections per star. The number of planet detections 
per time however is a function of target brightness. 
In practice, the effects of duty-cycle efficiency with 
respect to instrument calibration may further en- 
hance the appeal of massive stars. 

(4) Our simulations suggest that imaging planets 
that resemble the RV population will be unlikely 
around low-mass stars. We have considered a case 
where stars were chosen based on their age and 
brightness. However, another option to further 
improve the likelihood for d iscovery is t o conduct 
a metallicity-biased survey. lAgoll (|2007l ) discusses 
the prospects for such a strategy using space- 
based coronagraphs. Likewise, searches based on 
the exist ence of long-term RV tren ds may also be 
efficient (jCrepp and Johnsonll201ll ). 



(5) We have applied number statistics for the age 
of FG stars to K-type and early M-type stars. The 
discovery of significantly more young, low-mass 
stars would improve the prospects for imaging 
planets in the immediate solar neighborhood. Ob- 
servers are currently using activity-based youth 
indic ators to identify more nearby low-m ass tar- 
gets (|Liu et al.ll20"oalShkolnik et alj|2009h . 



(6) We have extended the iJohnson et al.l (|2010al ) 
results to stellar masses beyond M* = 2.OM , 
the range considered in their study. The planet 
occurrence rate may, however, saturate at high 
stellar mass. In this case, our calculations of the 
relative dctcctability for a = 1.0 and /3 = 1.2 
would over-estimate the relative promise of the 
most massive A-stars. However, given the number 
of super- Jupiters (m p sin i > 5Mj) discovered 
around massive K-giants with M* > 2.0M Q , it 
do es not appear t hat th e relationship measured 
bv IJohnson et al. ( 2010a) saturates i n the range 
2.0 < M*/M Q < 3 (|Sato et al"1l2010l ). 



but t his has yet to be measured (jRaghavan et al 
2010h . The number of planet detections depends 



(7) One challenge with measuring and interpreting 
planet detection rates is that dynamical effects, 
such as the gravitational impulse from nearby 
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passing stars, can deplete wide separation planets 
at young ages (jVeras et al Such encounters 

would preferentially effect l ow mass stars based o n 
binding energy arguments (jWeinberg et al.lll987l) . 
In this case, detection rates around MK stars may 
be artificially low and mimic the results for an al- 
ternative formation scenario. 

5. Summary and Concluding Remarks 

We have used Monte Carlo simulations to es- 
timate the number of extrasolar planets that are 
directly detectable in the solar-neighborhood as a 
function of stellar mass with current and next- 
generation high-contrast imaging instruments. 
Starting with the most recent empirical findings 
for gas-giant planet occurrence-rates from Doppler 
RV surveys, we consider ground-based observa- 
tions that survey volume-limited, and age-and- 
brightness selected samples of synthetic stars. Our 
calculations show that the probability for imaging 
a planet in the solar neighborhood lies heavily in 
favor of A-type stars as a result of strong corre- 
lations between star mass and planet properties. 
Such correlations are known to be present in the 
RV sample, and are required to reproduce direct 
imaging discoveries and non-discoveries to date 
(Fig. 2-5, 7). This result is independent of planet 
evolutionary model choice. It also holds when 
targets are pre-selected for youth and brightness, 
indicating that massive stars are ideal targets for 
reasons in addition to their natural youth and 
ability to serve as a bright beacon for AO correc- 
tion. 

The same effects responsible for creating a mul- 
titude of detectable planets around massive stars 
conspire to reduce the number orbiting low-mass 
stars. Loosening target selection criterion to in- 
corporate older and fainter stars, such as in the 
case of a volume-limited survey, can increase the 
number of planet detections around MK stars by a 
factor of several but at the expense of significantly 
lower efficiency. An AO system with a wave-front 
sensor operating in either of the I,Y,J bands can 
boost the number and rate of detections in this 
regime, suffering only a minor penalty with more 
massive stars. 

Using simulation input parameters that best re- 
produce observations in the solar neighborhood, 
we predict that surveys using instruments that 



generate contrast levels ofCo~5xlO~ 6 will pro- 
vide only a marginal increase to the known planet 
population at large separations, at most doubling 
the current number of detections. Surveys us- 
ing instruments with "extreme" AO that generate 
contrast levels of Co ~ 5 x 10~ 7 will provide a sub- 
stantial gain, permitting statistical analyses with 
a sample of 12-81 new planet imaging discoveries 
(Fig. 8). 

To maximize the planet yield, we find it effi- 
cient to observe massive stars at distances beyond 
50 pc, the canonical limit at which many observing 
campaigns are currently designed. Assuming non- 
ideal conditions, such as a bright sky background 
and decreasing planet semimajor axis distribution 
that truncates at 35 AU, dedicated observations 
can provide another factor of 2-3 in the number of 
directly imaged planets by extending survey selec- 
tion to greater distances, for example detecting at 
least one companion for every 17 stars observed 
with mass, 2.0 < M*/M < 2.5, out to 87 pc 
with C = 5 x 10~ 7 (Fig. 9). Searching closer to 
target stars enhances the number of detections in 
this regime by ks 11% when going from an IWA of 
4 X/D to 3 X/D. 

Stellar clusters represent a special case in high- 
contrast imaging because the constituent members 
have the same age and distance, neutralizing two 
important trade-offs between high-mass stars and 
low-mass stars. In this case, the number of de- 
tections is governed primarily by the stellar IMF, 
favoring low-mass stars by default (Fig. 10). If no 
correlations between star mass and planet prop- 
erties exist, the detection rate becomes flat as a 
function of stellar mass, and may have a small neg- 
ative slope because low-mass stars provide access 
to lower-mass planets. Simulations of a collection 
of stars resembling Tucana-Horologium show that 
this is the only case where M-stars can provide the 
most detections and also have the highest detec- 
tion efficiencies. The inclusion of any correlations 
whatsoever between star mass and planet proper- 
ties creates a positive slope, shifting the likelihood 
for detection back in favor of high-mass stars. 

At the same time, most nearby stellar clusters 
suffer from small number statistics, making it diffi- 
cult to image a planet around its high-mass mem- 
bers, because the stellar IMF falls off faster than 
the planet occurrence rate grows. This effect is ex- 
acerbated by the fact that: half of cluster members 
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are binaries and high-contrast programs usually 
avoid binaries; most stellar clusters are distant, 
«140 pc from the Sun, forcing the overall planet 
occurrence rate to take on its value at large separa- 
tions; and the occurrence rate of ultra-wide sepa- 
ration companions appears to be very small, of or- 
der several percent. Thus, MK stars will continue 
to dominate the number of planet detections in 
distant stellar clusters containing <200 members 
irrespective of formation mechanism^) Since ob- 
servers have generally either conducted age-and- 
brightness selected surveys restricted to the so- 
lar neighborhood, or systematic observations of 
nearby clusters, the current direct imaging planet 
population forms a dichotomy between early and 
late- type stars (Table 4). We find that this di- 
chotomy can be parsed according to the ratio of 
the effective spatial size of the young stellar group 
to the average distance of its members from the 
Sun. 

The planets orbiting HR 8799 and (3 Pic have 
much smaller semimajor axes and masses than 
those found in stellar clusters, and more closely 
resemble the RV population. Further suggestive 
of similarity is the fact that planets have been 
found around these two A-stars, while there ex- 
ists a dearth of detections around FGKM stars in 
the solar neighborhood despite dedicated efforts. 
The higher detection rate around A-stars com- 
pared to less massive stars is consistent with the 
correlations between star mass and planet proper- 
ties found in the RV population. 

Since the detection of planets orbiting HR 8799 
and /3 Pictoris requires that strong correlations ex- 
ist between star and planet properties, we propose 
that observers measure the slope of the detection 
rate curve for MK star cluster members. A flat 
or negative slope would indicate a separate forma- 
tion channel, specifically one with a weak corre- 
spondence between star and planet properties, ex- 
plaining the planets orbiting AB Pic, 1RSXJ1609, 
GSC 06214 and others. Such a diagnostic provides 
a complementary test to those involving measure- 
ment of metallicity correlations, semimajor axis 
distributions, and absolute occurrence rates. 



Although planets will eventually be imaged 
around nearby FG stars, under no circumstances 
do we find that intermediate spectral-type hosts 
have the highest detection rates. The results are 
either skewed to one side or the other. Along with 
the aforementioned biases, this may help account 
for the current list of exoplanet imaging discover- 
ies. Observations of A and M-stars can therefore 
provide the most leverage for discriminating be- 
tween planet formation scenarios. Surveys are cur- 
rently planned to specifically target stars in these 
mass regimes, using "extreme" AO at Palomar 
( Hinklev et al.ll2011bl ) and deep AO observations 
at Keck (Bowler 2011, private communication). 

In closing, we note that the first handful of 
Doppler-detected planets showed a strong prefer- 
ence for format ion around star s with an elevated 
metal content (jGonzalea I1997T ) . Subsequent ra- 



9 Notice however that aperture masking observations can 
access small star-planet separations and help to circum- 
vent this bias by targeting the bright member s of clusters 
l lKraus and I reland 2011; Hinklcy ct al. 201 laT). 



dial velocity observations have filled in this param- 
eter space, quantifying the relationship between 
planet occurrence and host star metallicity using a 
statistically significant sample dSantos et al. 2003; 
Fischer and Valentil 120051 : Ijohnson et al.l l2010al ) 
If the first few direct planet detections are a hint 
for which parameters govern formation at wide 
separations, an analogous situation in the field of 
high-contrast imaging may transpire involving cor- 
relations with star mass. 
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